The multiple effects of gradient coupling on network synchronization 
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Recent studies have shown that synchronizability of complex networks can be significantly improved by 
asymmetric couplings, and increase of coupling gradient is always in favor of network synchronization. Here 
we argue and demonstrate that, for typical complex networks, there usually exists an optimal coupling gradient 
under which the maximum network synchronizability is achieved. After this optimal value, increase of coupling 
gradient could deteriorate synchronization. We attribute the suppression of network synchronization at large 
gradient to the phenomenon of network breaking, and find that, in comparing with sparsely connected homo- 
geneous networks, densely connected heterogeneous networks have the superiority of adopting large gradient. 
The findings are supported by indirect simulations of eigenvalue analysis and direct simulations of coupled 
nonidentical oscillator networks. 

PACS numbers: 05.45. Xt, 89.75.-k 
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Complex networks have attracted a great deal of interest 
since the discoveries of the small- world [ 1 ] and scale-free [01 
properties. Roughly, small-world networks are characterized 
by a locally highly regular connecting structure and a globally 
small network distance, while the defining characteristic of 
scale-free networks is a power-law distribution p (fc) ~ k^^ in 
the node degree. Signatures of small-world and scale-free net- 
works have been discovered in many natural and man-made 
systems L3i,4i, JJ, and they constitute the cornerstones of mod- 
ern network science. 

At a systems level, synchronization is one of the most com- 
mon dynamical processes. For instance, in biology, synchro- 
nization of oscillator networks is fundamental [6]. In a com- 
puter network designed for large scale, parallel computation, 
to achieve synchronous timing is essential. Recent studies 
of the synchronizability of complex networks have revealed 
that small-world and scale-free networks, due to their small 
network distances, are generally more synchronizable than 
regular networks |0, [Ml- A somewhat surprising finding is 
that a scale-free network, while having smaller network dis- 
tances than a small-world network of the same size, is actu- 
ally more difficult to synchronize [9]. Considering the ubiq- 
uity of scale-free networks and the importance of synchro- 
nization in network functions, the finding seems to have gen- 
erated a paradox. However, it is recently found that, with 
weighted and asymmetric couplings, the synchronizability of 
scale-free networks can be significantly improved and, in 
gen eral, can be much higher than the small-world networks 
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For a pair of connected nodes on the network, the mutual 
couplings between them are usually unbalanced. One direc- 
tion will weight over the other direction and generate a cou- 
pling gradient on the link. To enhance network synchroniza- 
tion, both the direction and the weight of coupling gradient 
should be properly set according to the network properties 



such as node degree |[Tol[Tlll and betweenness IT2I1 . In Refs. 
Hol [Till it has been shown that, by setting the coupling gra- 
dient flow from the higher-degree node to the smaller-degree 
node, the synchronizability of scale-free networks can be sig- 
nificantly improved and higher than that of homogeneous net- 
works. With the same scheme of gradient direction, in Refs. 
lfril[T6ll it has been shown that the synchronizability of scale- 
free networks can be further improved by increasing the gra- 
dient weight, and larger gradient in general assumes higher 
synchronizability. This enhancing role of coupling gradient 
is further highlighted in Ref. where nodes are proposed 
to be connected by only gradient couplings, i.e. the one-way- 
coupling configuration. 

While most of the previous studies are focusing on the en- 
hancing role of coupling gradient on synchronization, there 
are accumulating evidences showing the opposite: large gra- 
dient may also deteriorates synchronization. For example, in 
Ref. ifTol [Tsll it is observed that as gradient increases, net- 
work synchronizability is firstly enhanced, gradually reaching 
to its maximum at an optimal gradient; then, after this optimal 
value, increase of gradient will suppress synchronization. The 
similar phenomenon is also briefly reported in Ref. I Till , there 
it is found that the suppressing effect of large gradient has a 
close relation to the network parameters. Despite of these ob- 
servations, a detail study on the multiple effects of coupling 
gradient on network synchronization is still absent. 

In this paper, we will study in detail the multiple effects 
of coupling gradient on network synchronization, and investi- 
gate the problem of synchronization optimization in asymmet- 
rically coupled scale-free networks. Our main findings are the 
following. (1) The destructive role of large gradient comes 
from the phenomenon of network breaking, increasing gradi- 
ent will also increase the breaking probability. (2) In general, 
small and densely-connected heterogeneous networks have a 
lower breaking probability than large and sparsely-connected 
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FIG. 1 : A schematic plot on the transformation of a general symmet- 
rical network (a) to a weighted asymmetrical network (b) and to a 
directed tree-structure network (the degenerative case)(c) by increas- 
ing the gradient parameter (3 in Eq. (1) from to oo. The dashed 
lines in (c) denote the links which are gradually diminished due to 
the increased gradient effect. 



homogeneous networks. (3) While large gradient deteriorates 
the propensity of global synchronization, partial synchroniza- 
tion of node clusters is enhanced. 

We consider oscillator networks of the following form 
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where F(xi) governs the local dynamics of uncoupled node i, 
H(x) is the coupling function, e is the coupling strength, and 
Gji is an element of the coupling matrix G which takes form 
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with ki the degree of node i and A = {oij } the adjacency 
matrix of the network, j = 1 if nodes i and j are connected, 
aij = otherwise, and ai^i — 0. To keep the synchronization 
state a solution of the system, we choose Gi^i = 1. 

It is worthy to note that the parameter (3 in Eq. (|2]i mod- 
ulates both the direction and weight of the coupling gradient 
on each link lfl6ll . Statistically, if /3 > 0, gradient is flow- 
ing from the higher-degree node to the smaller-degree node; 
while if /3 > 0, gradient is flowing in the opposite direction. 
Meanwhile, by tuning the absolute value we are also able 
to control the gradient weight: larger \j3\ generates larger gra- 
dient 1 16]. Please also note that changing (3 does not change 
the total coupling cost of the network, it only redistribute the 
weight of the couplings. Besides the flexibility of gradient 
control, the coupling scheme of Eq. ^ is also representative 
to a variety of network models proposed in previous studies. 
For instance, the symmetric network model in Refs. [7, 9] can 
be realized by replacing G with A in Eq. ([T]i; the weighted 
asymmetric network model (constructed based on the infor- 
mation of node degree) in Ref. 1 10] can be realized by setting 
/3 = in Eq. d2]l; and the directed tree-structure network 
model in Ref. [13] in principle can be realized by setting 
/3 ^ oo in Eq. (|2]i. A schematic plot on the realization of 
these different network models by changing (3 is illustrated in 
Fig. 1. 



The limiting case of /3 ^ cxd in Eq. (|2|i is of special inter- 
est [Fig. 1(c)], since itrepresents the extreme situation of node 
connection: the one-way coupling format, i.e., Gi,j = — 1 and 
Gj,i — when kj > ki, and Gij — and Gj,i — —1 other- 
wise. In this case, each node only receives coupling from one 
of its neighbors who has the largest degree. (Strictly speak- 
ing, each node is receiving coupling from one of its neigh- 
bors who has the largest node-scalar hi — fcf J2iev- ' ^''■'^ 
Vi stands the set of neighboring nodes of i. However, for 
complex networks of null degree-correlation, statistically we 
have hi ki. See Ref. [16] for details on the definition 
of node scalar). For this one-way coupled network, if it is 
non-degenerative (connective), nodes will be organized into 
a unique tree-structure topology, with the largest-degree node 
of the network locates at the root. Under this coupling con- 
figuration, the network reaches its maximum synchronizabil- 
ity MM, 

manifesting the constructive role of large gradi- 
ent. However, if there are two or more large-degree nodes on 
the network which are not directly connected, like the case of 
Fig. 1, the network will be degenerative (disconnected) [Fig. 
1(c)]. In such a case, the network will break into several sub- 
networks of tree-structure, and at the root of each subnetwork 
locates a local largest-degree node. Once is broken, the net- 
work can never be globally synchronized whatever the cou- 
pling strength, manifesting the destructive role of large gradi- 
ent. 

To have a global picture on the gradient effects, we have 
investigated the variation of network synchronizability as a 
function of the gradient parameter (3 for different network 
topologies. The synchronizability of coupled network can be 
evaluated by the method of master stability function (MSF) 
lilt], which states that a network is generally more synchro- 
nizable when the spread of the eigenvalue spectrum of its 
coupling matrix is narrow. In particular, let = Ai < 
A2 . . . < Aat be the eigenvalue spectrum of the coupling ma- 
trix G. Then the smaller the ratio R = the more 
likely synchronous dynamics is to occur on the network. In 
general, the matrix constructed by Eq. ^ is asymmetric and 
its eigenvalues are complex. Noticing that the coupling ma- 
trix G can be written as G = QID^^, with L = D - A, 
D = diag{ki, k2, .../cw} the diagonal matrix of degrees, and 
Q = diag{l/J2jLijkj,...l/J2jLNjkj} the normaliza- 
tion factors on rows of G. From the following identity 



det{QLD^^ - XI) = dci{Q^/^DP/^LDf^/^Q^''^ - XI) (3) 

it is found that the eigenvalues of the asymmetric ma- 
trix G are equal to that of the symmetric matrix H = 
qi/2£)/3/22^£)/3/2qi/2^ which are real and nonnegative. 

To simulate, we generate scale-free network of = 2^" 
nodes and average degree (fc) = 6 by the generalized model 
introduced in Ref. Ill 711 . In this model, the degree exponent 
7 can be adjusted via a parameter B. Defining the new pref- 
erential attachment function as p ~ {ki + B)/J2j{kj + B), 
it can be proven that 7 — 3 + B/m, with m — 3 the num- 
ber of new links that associated to each new added node in 
the model. Using B = 0, 12 and 51, we generate scale-free 
networks of degree exponents 7 = 3, 7 and 20, respectively. 
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FIG. 2: (Color online) For scale-free networks of TV = 2 node 
and average degree (fc) = 6. (a) The variation of eigenratio R as 
a function of gradient parameter /3 for networks of degree exponent 
7 = 3 (the lower curve) and 7 = 7 (the middle curve) and 7 = 20 
(the upper curve). The optimal gradient /3o is about 0.9 in the later 
two cases. Each data is an average result of 100 network realizations, 
(b) For the case of 7 = 7 in (a), the value of Ri for different network 
realizations under parameters /3 = 0.5 (the red curve) and (3 = 5 
(the black curve). The extreme values of Ri indicate the broken of 
the network topology in the corresponding realizations. 



The variations of i? as a function of /3 for these three net- 
works are plotted in Fig. 2(a). It is found that, in the case 
of 7 = 3, increasing (3 will enhance synchronizability mono- 
tonically; however, in the case of 7 = 7 or 8, as f3 increases 
from zero, the network synchronizability is firstly enhanced, 
and, after reaching its maximum at value about /3o « 0.9 1 19], 
it begin to be suppressed. Another interesting finding in Fig. 
2(a) is that, for the fixed gradient parameter (3, increasing the 
degree exponent 7 will always decrease the network synchro- 
nizability, indicating the superior synchronizability of scale- 
free networks under gradient couplings. 

To gain insight on the transition of gradient effect from en- 
hancing to suppressing synchronization, we go on to investi- 
gate the changes happening in the neighboring region of the 
optimal gradient f3o. With 7 = 7 (the middle curve in Fig. 
2(a)), we plot in Fig. 2(b) the individual value of Ri for a 
large number of network realizations under gradient param- 
eters /?! = 0.5 < /3o and /Jg = 5 > /3o. It is found that, 
for /3i — 0.5, the eigenratio Bi is oscillating around its mean 
value i? w 6 with very small fluctuations; while for f3i = 5, 
the eigenratio Ri occasionally bursts into some extreme val- 
ues of order 10^. Since of i? = a divergent Ri thus 
indicates A2 0, which, from the eigenvalue analysis, im- 
plies the breaking of network topology in the corresponding 
realization. Therefore the suppression effect of large gradient 
can be attributed to the phenomenon of network breaking, and 
the optimal gradient can be understood as a balance between 
the enhancing and the suppressing effect of gradient coupling. 

By knowing that synchronization suppression at large gra- 
dient is induced by network breaking, we next to investigate 
the relationship between the breaking probability and the net- 
work parameters. To facilitate the analysis, we consider again 
the limiting case of /? ^ 00 in Eq. As shown in Fig. 
1(c), network breaking happens when there are more than 
one local-maximum-degree nodes coexist on the network. To 



FIG. 3: (Color online) For scale-free network of A'^ = 3 x 10^ nodes 
and under the limiting case of /3 00 in Eq. the probability of 
network breaking pt as a function of the degree exponent 7 (a) and 
the average degree (k) (b). Each data is an averaged result of 1000 
network realizations. 



break the network, the "breaking nodes" do not have to be 
possessing very large degree, they are only required to have 
the largest degree among its neighbors. Once is broken, the 
network will be divided into several tree-structured subnet- 
works, with each subnetwork is led by one of such "break- 
ing nodes". Due to the complicated configurations of the 
"breaking nodes", we are unable to give an analytical predic- 
tion on the relationship between the breaking probability and 
the network parameters. The numerical results on the vari- 
ation of pb as functions of 7 and (fc) are plotted in Fig. 3, 
which show that, in comparing sparsely connected homoge- 
neous networks, densely connected heterogeneous networks 
are more difficult to break down ifTsll . It can be also found 
from Fig. 3(b) that the increase of system size will always 
increase the breaking probability. These findings (calculated 
for the limiting case) are coincident with the findings in Fig. 
2 (calculated for the general case), both tell us that heteroge- 
neous networks are more sustainable to large gradient. 

The results exemplified in Figs. 2 and 3 are based on 
the eigenvalue and topology analysis. It is useful to exam- 
ine the gradient effects in real oscillator networks. For this 
purpose we have check the synchronization of scale-free net- 
works of coupled nonidentical chaotic Rossler oscillators, a 
typical model employed in studying network synchronization 
iS Hi m [H. The dynamics of a single oscillator is de- 
scribed by Fj(xj) = [-ojiyi - Zi,uj.iXi + 0.15yi,Zi(xi - 
8.5) + 0.4], where uji is the natural frequency of the ith oscil- 
lator. In simulations we choose uJi randomly from the range 
[0.9, 1.1], so as to make the oscillators nonidentical. The cou- 
pling function is chosen to be H(x) = x. The degree of 
synchronization can be characterized by monitoring the am- 
plitude A of the mean field X{t) = '^fLi x.Xt)/N |10]. For 
small coupling strength e, X{t) oscillates irregularly and A 
is approximately zero, indicating lack or a lower degree of 
synchronization. As the coupling parameter is increased, syn- 
chronization sets in and A is increased gradually from zero 
(nonsynchronous state) to its maximum (synchronous state). 
By e = 0.15, we plot in Fig. 4(a) the variation of A as a func- 
tion of (5. It is found that as /? increases A is firstly increased, 
and reaching its maximum at about fio ~ 1, manifesting the 
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FIG. 4: For scale-free network of = 2 x 10 nonidentical R6 
ssler oscillators and average degree {k) = 6. (a) For e = 0.15, the 
variation of the mean field amplitude A as a function of the gradient 
parameter (3. Network synchronization is optimized at about /3 ~ 1. 
Each data is an averaged result of 50 network realizations. For a 
chosen network realization of (a), the frequency distribution of the 
oscillators under gradient parameters /3 = 0.1 (b), /3 = 1 (c) and 
/3 = 10 (d). 



constructive role coupling gradient. Then, as /3 increases from 
Po, A begins to decrease, manifesting the destructive role of 
coupling gradient. An interesting finding is that increasing f3 
further does not decreases A continuously, the system always 
keeps on high coherence at about A « 8. Please note that 
for the adopted coupling function H, the stable region of the 
MSF function |8] has only a lower boundary, which confirms 
that the decrease of A at f3 > f3o in Fig. 4(b) is induced exclu- 
sively by the breaking effect (not induced by the instability of 
the shortest wave mode |8]). 

The finding that A keeps on large values at very large (3 
indicates that, despite of the increased probability of network 



breaking, nodes are still strongly correlated, with a manner 
that is different to the situation of small j3. To gain insight, 
we choose a network realization in Fig. 4(a), and plot the fre- 
quency distribution uji of the oscillators under different gra- 
dient parameters: (3 = 0.1 in Fig. 4(b), f3 ~ 1 in Fig. 4(c) 
and (3 ~ 10 in Fig. 4(d). It is found that, for small gradi- 
ent [Fig. 4(b)], LUi is distributed randomly around the mean 
value u! — 1, indicating the low system coherence; around the 
optimal gradient [Fig. 4(c)], most of the oscillators are syn- 
chronized to be having the same frequency, with few excep- 
tions which are synchronized to another frequency, indicating 
a higher coherence of the system dynamics under this gradi- 
ent; for large gradient [Fig. 4(d)], the oscillators are separated 
into two clusters of the similar size, the frequency of one clus- 
ter is different to that of the other one. Fig. 4(d) indicates that, 
under the large gradient, although the probability of network 
breaking is high, network nodes are still strongly correlated 
due to the existence of synchronous clusters. The clusters are 
a direct result of the network breaking. Led by the "breaking 
nodes", each subnetwork develops into a synchronous cluster. 

In summary, we have studied the multiple effects of cou- 
pling gradient on network synchronization, and investigated 
the dependence of these effects to the network parameters. 
Our findings suggest that, in comparing with sparsely con- 
nected homogeneous networks, densely connected heteroge- 
neous networks take more advantages from the gradient cou- 
plings. 
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